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We present the calculation of heavy Wilson quark mass renormalization constants from large beta Monte Carlo 
simulations. Simulations were performed at various beta larger than 9, each on several spatial lattice sizes to 
allow for an infinite volume extrapolation. We use twisted boundary conditions to suppress tunneling and work 
in Coulomb gauge with appropriate adjustments for the temporal links. The one-loop coefficient obtained from 
this method is in agreement with the analytical result and a preliminary result for the second order coefficient is 
reported. 



1. INTRODUCTION 

As actions for lattice calculations become more 
complicated, short distance calculations for nor- 
malizing the operators in the action become more 
and more difficult. It is useful to have as many 
tools as possible to attack the challenge. Several 
types of automated calculation of Feynman dia- 
grams and nonperturbative methods have been 
developed. Another alternative is the use of di- 
rect calculation of short distance quark-giuon 
Green functions by Monte Carlo methods . Last 
year this method was applied to a determination 
of the additive mass renormalization for static 
quarks to third order ||]. In this paper, we present 
a preliminary estimate of the mass renormaliza- 
tion for heavy Wilson quarks. 

2. SIMULATION 

There are several potential problems in naively 
performing simulations at large f3 to extract per- 
turbative information. One is the problem as- 
sociated with infra-red divergences coming from 
the zero modes. The second is the phenomenon 
of tunnehng between the various vacua of SU{3). 
The third is finite volume effects, since the phys- 
ical sizes of space-time used in the simulations 
are quite small. All of these effects may be re- 
moved or reduced by using twisted boundary con- 
ditionsQ for the fields. These special boundary 
conditions remove the zero modes completely and 



suppresses tunneling for sufhciently large 
The remaining finite size effects can be removed 
by extrapolating to the infinite volume limit. 

2.1. Twisted Boundary Conditions 

We use a 3-way periodic twist (spatial direc- 
tions only) for the gauge fields defined by. 
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where 77 = e^'^*/^ and = ^y^l- The rj's are 
constant SU{3) matrices. 

2.2. Gauge Fixing 

We must fix the gauge in order to extract en- 
ergies from quark propagators. Our gauge fixing 
procedure is as follows: 

(1) fix to Coulomb gauge, then 

(2) fix the temporal links by maximizing, 
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where Z(t) are elements of the center of S'C/(3). 

2.3. Parameters 

The following values of the coupling were cho- 
sen: /3 = 9, 10, 12, 15, 20, 40 and 60. Lattice sizes 
of Ls = {4, 6, 8, 10, 12, 14, 16} and U = 16 were 
used at each coupling to allow for an extrapo- 
lation to the infinite volume limit. In this pre- 
liminary study, we have accumulated 500 2000 



configurations separated by 100 sweeps. The er- 
rors after binning did not show significant residual 
correlations. 

3. STATIC QUARK SELF-ENERGY 

We test our gauge-fixed configurations by cal- 
culating the first two coeSicients of the static 
quark self-energy and compare to analytic|3] 
and previous high /3 Monte-Carlo calculations] 2] 
where no gauge fixing was performed. We also 
compare results from the Polyakov loop and the 
static quark propagator. 

3.1. Polyakov Loop 

The tadpole-improved static quark self-energy 
(in a finite volume) was calculated from the mea- 
sured value of the Polyakov loop according to, 
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where uq is the mean link from the average pla- 
quette. Our fitting procedure to extract the 
renormalization constants is as follows: 

1. Fit the measured values at fixed i to a poly- 
nomial in av{q*) to obtain ci{L) and C2{L), 
i.e. Ea{L)/a = ci(i) + ac2{L) + a^C3(L). 

2. Extrapolate the coefficient ci{L) using 
ci{L) = ci + Xi/L to get the one- loop co- 
efficient, Ci. 

3. Extrapolate C2(L) with C2(£) = C2+X2/ 
l3oXiln{L)/L where Po = ll/27r to get the 
two loop coefficient, C2. 

The functional form and the coefficient of the 
logarithm comes from the anticipation of the 
(running) Coulomb interaction with the image 
charges with some effective distance L/Xi. 

We also perform the fit in the reverse order, 
i.e. extrapolate to infinite volume at each f3 and 
then fit to a cubic polynomial in ay. In the 
infinite volume extrapolation, we use Eq{L) ~ 
Eo + Yi/L + (3oXiln{L)/L where Xi is the co- 
efficient from the previous fit. The results were 
in perfect agreement. The renormalized cou- 
pling in the V-scheme[0, ay, was evaluated at 
q* = 0.84a-i. 





Cl 


C2 


analytic 
static propagator 
Polyakov loop 
Reference 


1.070 
1.071(1) 
1.070(3) 
1.070(3) 


0.11 
0.20(7) 
0.17(14) 
0.14(4) 



Table 1 

Results for the static quark self-energy. 

3.2. Static Quark Propagator 

On the same set of lattices, we have measured 
the static quark propagator averaged over the 
spatial sites. These were then fit to single ex- 
ponentials starting from timeslice 2 for which we 
get a good fit. The perturbative coefficients were 
extracted in the exact same manner as in the pre- 
vious section. 

3.3. Comparison 

We compare our results with the analytical cal- 
culations of Heller and KarschQ and with previ- 
ous measurements Q of the first two coefficients in 
Table |l|. Our measurements on gauge fixed con- 
figurations with both gauge invariant and covari- 
ant operators are in good agreement with those 
in the literature. These results give us confidence 
in our gauge fixing procedure. 

4. WILSON QUARK REST MASS 

We now apply similar methods to the Wilson 
quark mass renormalization. In order to accom- 
modate twisted boundary conditions consistently 
for propagating quarks, we enlarge the quark in- 
ternal degrees of freedom to Nspin x Ncoior x 
Nsrneii whcrc the smell group is SU{Nsrneii = 
3) Q . The boundary condition for the quark field 
reads, 

'f{x + LjCj) = %*(x)17Tie*'^/^-'<"- (4) 

where the phase factor ensures anti-periodicity. 
The momenta of the quarks are then given by, 

P. = Tr^^(2n. + i) (5) 



* coior^i 



- Uy + = mNcoior (6) 

where m and Ui are integers. The form of the 
constraint on is due to our particular choice of 



the twist matrix in the z-direction. 

We have measured the ground state only, (rix — 
0, Uy = 0, = 0), for which the color- smell wave- 
function is proportional to the identity matrix 
The extrapolation to the infinite volume limit will 
then give the rest mass, E{p — 0), 



where Afj"' = + amo), aniQ 
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1/(2k)-4 and 
A/j^' is calculated in Ref. ||]. As explained in 
Ref. 1^, the parameters in the Wilson and clover 
actions should be taken as mass-dependent when 
the masses are large. Here we calculate the mass 
renormalization for k = 0.1. 

The ground state momenta was projected out 
both at the source and the sink. Only the spin 
components and 2 (set equal to each other) were 
used in the simulation. 

4.1. Fitting Procedure 

We have fitted the quark propagator to a cosh 
starting at time-slice 3 for which we get a good 
X^. The tree- level mass, Afj°'(L), was then sub- 
tracted from the raw data so as to cancel some 
of the finite volume dependence of the rest mass. 
The extraction of the one-loop coefficient, Afj^', 

[21 

and the two-loop coefhcient, Afj , was done using 
the same procedure as the static quarks. How- 
ever, we have found that a and XjL^ term 
was required to get a good fit to the one-loop co- 
efficient (see Fig. H). The second order coefficient 
was extracted with a 1/L and ln{L)/L term. 

As in the case of the static energy, we fit the 
masses in the reverse order to check for consis- 
tency. The infinite volume extrapolation is done 
with 1/L, l/i^^ 1 /i3^ ^j^]^ ^-j^g coefficient 

of ln{L)/L fixed as before. We show one example 
in Fig. ||. The fit to a third order polynomial in 
ay is shown in Fig. ^ where the value of q* was 
obtained in Ref. . The results are tabulated in 
Table |. 

5. SUMMARY 

It was shown that the one-loop Wilson quark 
mass renormalization constant is reproducable 
with a combined extrapolation and statistical er- 
ror of 0.5% from large (3 Monte-Carlo simulations 
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Figure 1. The infinite volume extrapolation of 
Ml - m[°^ at P = 60. Note that the point L = 24 
was not used in the fit. 
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□ M['l = 3.8275 analytical 
m('1 = 3. 827(18) fit 
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Figure 2. The one- loop coefficient as a function 
of the spatial lattice size. 
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Figure 3. The two-loop coefficient as a function 
of the spatial lattice size. 



of the quark propagator. Twisted boundary con- 
ditions ensured a perturbative expansion around 
the trivial vacuum and the absence of zero modes. 
The gauge fixing procedure was checked to repro- 
duce the static quark self-energy using both the 
Polyakov loop and the static quark propagator. 

Our preliminary result for the two-loop Wil- 
son quark mass renormalization (for k — 0.1) is 
—7(1) (0.5). We have a crude result for the three- 
loop term (see Fig. ^), but it is poorly determined 
and sensitive to the details of the fit procedure. 
For now, we simply take it as an estimate of the 
systematic error in the determination of the two- 
loop term. A more detailed study of the sys- 
tematic errors is now in progress including the 
infinite volume extrapolation, dependence on q* 
and truncation of the series. The heavy quark 
mass renormalization for improved actions, such 
as clover and NRQCD, are also currently under 
investigation. 
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□ 3.8275 : B.Mertens et al 
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Figure 4. The fit to {Mi-M{')/av of the infinite 
volume extrapolated masses, {q* = 1.75a~^) 









analytic 
mW(L) extrapolation 
reverse order 


3.8275 
3.827(18) 
3.827(18) 


unknown 
-7.6(1) 
-7(1) 


our estimate 


3.83(2) 


-7(1) 



Table 2 

Results of the fit to the Wilson quark mass renor- 
malization constants for k = 0.10 
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